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Sunmary 

Using  the  fact  that  cavities  are,  at  low  cavitation  numbers,  approximately 
prolate  spheroids,  this  report  develops  the  theory  of  wall  corrections  in  closed 
tunnels  and  free  jets  of  circular  cross-section. 

The  results  are  displayed  graphically  for  a wide  range  of  the  ratio  ( cavity 
/ length)  / ( tunnel  width) , 

In  addition,  a first-order  theory  is  developed  which  applies  only  to  small 
cavities  at  low  cavitation  numbers.  Hie  wall  correction  is  then  shown  to  be 
proportional  to  the  cube  of  this  ratio,  whereas  in  two-dimensional  theory  it  is 
proportional  to  the  square  of  the  ratio.  Moreover , the  corrections  in  a dosed 
tunnel  an  shown  to  be  roughly  four  times  as  great  as,  and  in  the  opposite  sense 
to,  those  applicable  to  free  Jets.  In  two-dimensional  theory  they  are  twice  as 
great. 

The  unaided  theory  yields  information  on  the  thickness  ratio  only  of  the 
cavity.  However,  if  the  additional  assumption  is  mads  that  the  drag  coefficient 
is  relatively  insensitive  to  boundary  effects  - an  assumption  which  is  oertainly 
true  for  the  two-dimensional  oaae  - than  the  length  and  breadth  of  the  cavity 
may  be  investigated  separately.  It  is  thus  found  that,  in  the  axi symmetric  as 
in  the  two-dimensional  oase,  the  boundary  correction  to  the  length  of  a small 
cavity  is  proportionately  twice  as  great  as  that  to  the  width. 
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Symbol* 

Homan  letters 


A, B  auxiliary  functions  defined  in  equations  (7.1). 

* 

o the  radius  or  semi -width  of  the  cavitating  head. 

Ct ( x ' ) an  auxiliary  funotion  of  used  in  equation  (4.14). 

Cj_)  the  true  drag  coefficient,  based  on  the  area  of  aross-section  of  the 

o*  vita  ting  head. 

Oq  the  modified  drag  coefficient,  based  on  the  area  of  cross-section  of 

the  cavity. 

C„  „ coefficients  defined  in  Appendix  II. 

0 *(*')  an  auxiliary  funotion  of  r' used  in  equation  (4.14). 

B,  %'  incomplete  elliptic  integrals  defined  in  equation  (4.5). 

f the  semi-length  of  the  distribution  of  sources  and  sinks  on  the  axis. 

P,  P'  inoonplete  elliptic  integrals  defined  in  equation  (4.5). 

' g(  6)  an  auxiliary  function  defined  in  equation  (4.2). 

C a function  relating  the  cavitation  number  to  the  sise  of  the  cavity. 

I(jio)  an  auxiliary  function  defined  in  equation  (6.3). 

J v the  Bessel  function  of  the  first  kind  and  of  order  v . 

k,  k'  the  modulus  and  co-aodulus  of  elliptic  integrals  defined  in  equation 

(4.4).  Physically,  k'  provides  a rough  measure  of  the  cavitation 
number. 

ky^y  the  m th.  positive  zero  of  Jy(z),  taken  in  ascending  numerical  order. 

K a cooiplete  elliptic  integral  defined  in  equation  (4.5). 

a nU  is  the  axial  component  of  fluid  velocity. 

m«  the  particular  value  of  m at  the  equator  of  the  cavity, 

M a function  relating  the  strength  of  the  uniform  stream  to  the  size  of 

the  oavlty. 

p a parameter  occurring  in  equation  (3.2). 

Pq  the  uniform  pressure  along  the  cavity  wall. 

p„  the  pressure  in  the  undisturbed  etream,  far  from  the  oavlty. 

q a mathematical  parameter  introduced  in  equations  (4.14).  Physically, 

it  providss  a rough  measure  of  th*  ratio  of  the  size  of  the  cavity 
to  the  diameter  of  the  etream. 

Q cavitation  number. 

f 

I 
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■ G ( u m ) , tha  oavitation  number  in  tha  unbound  ad  stream  which 
would  yield  a cavity  of  tha  length  actually  observed  in  tha  bounded 
stream. 


« G ( u o ) , tha  oavltation  number  in  tha  unbounded  stream,  which 
would  yield  a cavity  of  tha  width  aotually  observed  in  the  bounded 
stream. 

denotes  tt  remainder  term  of  smaller  order  than  the  terms  retained, 
a variable  of  summation, 
a variable  of  integration, 
velocity  of  the  undisturbed  stream. 

non-dimensional  coordinate  parallel  to  axis  of  eynnetry)  unit  of 

) length  is 
) radius  of 

non-dimensional  coordinate  perpendicular  to  axis  of  ) tunnel 
symmetry  ) 

the  semi-length  of  the  cavity. 

the  equatorial  radius  or  semi-width  of  the  cavity. 


auxiliary  functions  of  *r'  used  in  equations  (4.14). 


Ireek  Letters 

r, , T, , r,  auxiliary  functions  of  q,  defined  in  equations  (4.17). 

5 a mathematical  parameter  introduoed  by  equations  (4.1).  Ehysioally, 

i t provides  a rough  measure  of  the  ratio  of  the  else  of  the  cavity 
to  the  diameter  of  the  stream. 

t an  operator  i ml  eating  the  difference  between  a bounded  stream  value 

and  the  corresponding  unbounded  stream  value. 

A i , A* , 4, , d4  terms  contributing  to  A I ( g 0 ) • 

e the  displacement  of  the  jet  surface  at  the  plane  of  eyonetry. 

k‘  a mathematical  parameter  defined  in  equation  (4.15)*  Physically,  it 

is  a function  of  the  cavitation  number  only. 

\fu  ccniooidal  coordinates  defined  by  equation  (1.1). 

u 0f  um  particular  values  of  g at  the  equator  and  pole  of  the  oavity, 

respectively. 

v a parameter  taking  the  value  0 far  a free  surface  and  1 far  a fixed 

boundary. 

p the  x -coordinate  of  a source  element. 


u 


V*'  )’> 
Yt(  )J 


p density  of  ttM  liquid, 

^(x,  y)  that  port  of  the  axial  conponont  of  velocity  oontributed  by  the  lino 
aouroe-aink. 

, d (x^.w.p)  * potontlol  function  dofinod  in  equation  (3*9)* 
f(x,y)  Stoke*'  a trow  function. 

Q(x,y)  a current  function  dofinod  in  aquation  (1.5). 

(fc(x,y)  that  part  of  the  current  function  oontributed  by  the  lino  acuroe-aink. 

0(x,y,i/,p)  a particular  current  funotion  defined  in  aquation  (3.2). 


A bar  indloatea  that  the  ayafcol*  apply  to  the  unbounded  stream. 


Introduction 


Tho  purpoao  of  this  re port  la  to  investigate  mil  offaota  on  the  overall 
proportion*  of  a oavity  formed  by  an  axially  aye  trio  body  in  rapid  notion 
relative  to  a stream  of  eater.  In  the  laboratory  it  ia  poaaible  to  naaaure  the 
dinenaiona  of  auoh  a cavity , formed  in  a stream  having  * certain  finite  oroaa- 
sectional  area.  The  problem  then  ia  to  calculate  what  the  dimsnaiona  would  be 
if  such  a oavity  were  formed  in  a atreaa  of  unlimited  extent  in  all  directions. 

There  already  exiata  a oonaid arable  literature  on  the  theory  of  mil  oorreo- 
tiona  applicable  to  non -oa  vita  ting  flow  in  wind-  or  mter-tunnela.  In  that  case, 
however,  the  boundariea  of  the  flow  are  determined,  being  the  walla  of  the  tunnel 
and  the  surface  of  the  aerodynamio  body  under  test,  and  the  problem  consists 
of  determining  the  oorreotiona  to  be  applied  to  the  veloci ty  or  pressure  field. 

In  our  oase,  on  the  other  hand,  the  internal  boundary  of  the  flow  is  variable, 
being  the  surfaoe  of  the  oavity,  and  only  the  pressure  distribution  aloi^  it 
is  known. 

In  two  dimsnaiona,  tbs  problem  has  been  solved  for  oertain  simple  head 
shapes  by  the  use  of  oonforaal  transformations  (ref.  1 aid  6).  In  three 
dimensions,  even  under  the  assumption  of  axial  synsaotry,  this  powerful  method 
is  no  longer  available.  ▲ partial  and  approximate  solution  of  the  svisliy 
symmetric  problem  will,  however,  be  derived  in  the  following  pages  by  the  use 
of  aouroe-sink  methods. 

The  basic  fact  underlying  our  present  method  is  that  the  shapes  of  such 
cavities  are  known  to  apprniflmdhi  to  prolate  spheroids.  Furthermore , it  ia 
poaaible  to  produce  closed  stream  surfaoe  a,  whioh  are  either  exactly  or 
approximately  spheroidal  in  shape,  by  combining  a very  sliqple  axial  aouroe- 
sink  distribution  with  a uniform  stream  of  either  unbounded  or  bounded  extent, 
respectively. 

The  characteristic  feature  of  oavity  flow  ia,  of  course,  that  the  fluid 
pressure  is  uniform  along  the  free  surface  of  the  oavity,  and  a dimensionless 
oavitation  number,  Q,  is  accordingly  defined  by  the  ratio 


<J  - 


P«“Po 


(0.1) 


Here  p^  and  pQ  are  pressure  in  the  indlsturbed  stream  and  the  pressure  on 
the  oavity  wall  respectively,  whilst  p and  17  are  the  density  of  the  liquid 
and  the  velocity  of  the  undisturbed  stream. 

In  the  flow-patterns  which  we  shall  construct  by  source-sink  methods,  the 
pressure  along  the  closed  stream  surfaoe  will  be  only  approximately  uniform, 
this  constituting  an  intrinsic  imperfection  of  the  method,  and  we  shall  adopt 
the  oonvention  of  considering  the  pressure  at  the  equator  of  the  stream- 
surfaoe  as  characteristic  of  the  "oavity".  The  oavitation  mmfcer  will  then 
be  defined  In  terms  of  this  equatorial  pressure. 

Without  ary  further  assumptions  we  can  then  investigate  the  effect  of  the 
boundary  on  the  thickness  ratio  of  the  oavity,  at  a given  oavitation  number. 

In  order  to  invest!^ te  the  behaviour  of  the  length  and  width  separately, 
however,  we  mrta>  the  as*usg>tion,  rendered  plausible  by  the  known  results  of 
two-dimensional  theory,  that  the  drag  ooaffioient  of  a small  oavity  is 
unaffected,  to  a high  degree  of  accuracy,  by  the  boundary. 
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1.  Notation 


In  tills  report  it  will  be  convenient  to  employ  two  non-dimensional 
coordinate  systems,  an£  the  necessary  transformations  are  set  out  hare  for 
ease  of  reference.  The  cylindrical  coordinate  system  (Pig.  l)  will  employ 

x,  measured  along  the  axis  of  qrasetry  from  the  centre  of  the  oavity,  and 

y,  measured  perpendicularly  from  the  axis  of  aymsetry.  The  oonlooidal 
coordinate  system  (Pig.  2)  will  enplay  X , which  is  oonstant  on  each  member 
of  a oonfoodl  system  of  hyperboloids  of  revolution,  and  n , which  is  oonstant 
on  each  member  of  the  orthogonal  gystem  of  oonfooal  prolate  spheroids.  The 
fool  of  this  aye  tea  are  at  y ■ 0.  x ■ * f , where  f is  a positive  parameter. 
The  relations  between  (x,y)  and  ( x,/i)— are 


x/t  ■ tanh  X ooth  u 
y/f  a aeoh  X ooseoh  u ; 


(1.1) 


and,  in  particular, 

alien  x ■ 0 X » 0 , 

ldian  y a 0 (|x|  > f ) , 

Aon  y a 0 (|x|  < f)  u « ••  . 

We  ahall  study  the  shape  of  the  closed  stream-surface  generated  by  the 
superposition  of  a uniform  stream  of  velocity  U,  parallel  to  the  x-axla,on  a 
certain  source -oink  distribution.  The  latter  extends  along  the  x-axis  from 
- f to  ♦ f and  has  a density  equal  to  x.  The  unit  of  length  is  the  radius  of 
the  stream  at  infinity,  so  that  the  limiting  oaae  of  an  unbounded  stream  is 
obtained  by  letting  f tend  to  aero.  The  fineness  ratio  of  the  closed  stream- 
surface  is  controlled  by  the  choice  of  U. 


The  semi-length  of  the  closed  stream-surface,  that  is,  the  value  of  x 
when  y ■ 0,  will  be  denoted  by  xe.  Similarly,  the  equatorial  radius  of  the 
closed  stream  surface,  that  ia,  the  value  of  y when  x * 0,  will  be  denoted 
by  y0  . On  the  same  basis,  the  value  of  it  at  x ■ 0,  y ■ y0f  will  be  denoted 
by  , and  the  value  of  a at  x * x0,  y ■ 0,  will  be  denoted  by  um.  Thus 


y0  ft  • cosech  pe 
x0  /f  a coth 


(1.3) 


The  results  of  the  analysis  will  be  displayed  in  the  form  of  a comparison 
between  the  properties  of  the  dosed  stream  surface  for  general  values  of  f 
with  the  corresponding  properties  in  the  limiting  case  of  the  unbounded  stream. 
Aooardiigly,  it  will  be  convenient  to  use  a bar  to  indicate  values  applying  to 
the  oeee  of  the  unbounded  stream,  and  the  symbol  h to  indicate  the  difference 
between  the  value  for  the  bounded  stream  and  that  for  the  unbounded  streaau 
Thus,  for  exssple, 

0 ■ 0 ♦ dU  . (1.4) 


It  is  usual  to  describe  axiayometric  flow  patterns  by  means  of  Stokes' 
stream  function  *(x.y),  where  2w#(x,y)  is  the  rate  of  volume  flow  of  liquid 
through  the  oirole  (x,y)  in  the. positive  direction  of  x.  Because,  however, 
this  function  vanishes  identically  on  the  axis  of  eynnetry,  it  will  be  more 

convenient  here*  to  use  a current  function  (l(x,y)  defined  by 

• . 


y*  n(x,y)  « *(x,y)  . 


(1.5) 


This  function  hkf  ifL  obvious  physical  significance,  for.  it  is  equal  to  half 
the  average  x-4csq>unhnt  pi  telocity  inside  the  circle (x,y).»  It  vanishes  on 


the  x-eocls  onfcr^at ♦spwdi*T points , namaly  the  stagnation  patyt$  (•♦  *e.O).  v 

' ru  A,-:  -..t-V  * ■ -*■  ,v 


Tha  closed  stream  surface  «(x.y)  ■ 0 breaks  up  into  two  tranches, 
namely  the  axis  of  symmetry  y ■ 0 and  the  closed  stream  surfaos  0(x,y)  ■ 0. 

It  is  only  the  latter  branch  which  is  of  interest  in  this  paper. 

Prom  this  point  onwards,  in  the  interests  of  brevity,  we  shall  replace 
the  strictly  correct  appellation  "closed  stream  surfaoe"  by  the  rather 
loose  term  "cavity" . 

The  axial  oooponent  of  the  fluid  velooity  at  any  point  of  the  flow  will 
be  denoted  by  at),  so  that  m is  a lunotion  of  x and  y,  or  of  X and  u . In 
practice,  however,  we  shall  only  study  the  variation  of  m in  the  plane  x a 0, 
where,  by  symmetry,  mil  is  the  actual  magnitude  of  velocity.  The  particular 
value  of  m at  the  equator  of  the  cavity  will  be  denoted  by  m0,  and  it  follows 
easily  from  Bernouilli’s  law,  and  the  definition  (0.1)  of  the  cavitation 
number,  that 


1 ♦ Q ■ . (1.6) 

It  is  often  convenient  in  oavltaticnal  theory  to  consider  two  drag 
coefficients,  the  first,  Cq,  being  based  on  the  area  of  the  plate  which  Is 
supposed  to  be  causing  the  cavity;  and  the  second,  Cjj,  being  based  cn  the 
maxi  mm  cross-sectional  area  of  the  cavity.  It  will  be  convenient  to 
distinguish  between  the  two  by  referring  to  them  in  this  report  as  the  "true" 
and  the  "modified"  drag  coefficient,  respectively.  If  the  radius  of  the 
oavitating  plate  is  denoted  by  o,  then  we  have  the  obvious  relationship 

yo'/o*-  Cj/Cfl  . (1.7) 

Simmons  (ref.  1 ) has  shown,  from  considerations  of  momentum  and  continuity, 
that 

Crf  - Q - 2yZ*  I (m  - 1)*  y dy  , (1.8) 

■'y  0 

the  integral  being  evaluated  along  the  y axis  (x  ■ 0).  It  will  be  possible 
to  oalculate  the  boundary  effect  on  Oq,  and  hence,  by  making  a plausible 
assumption  as  to  the  behaviour  of  Cq,  it  will  be  possible  to  deduce  the 
boundary  effect  on  the  ratio  y^o  from  equation  (1.7). 

2.  Cavity  in  an  unbounded  stream 

ffe  develop  first  the  theory  for  the  unbounded  stream,  as  the  analysis 
required  is  simpler  than,  and  completely  dissimilar  from,  that  required  far 
the  stream  of  finite  lateral  extent.  The  fact  that  this  case  corresponds  to 
a vanishingly  w— 11  value  of  f is  no  inconvenienoe , since  all  the  ratios  of 
length  which  occur  in  this  section  are  finite  in  sise. 

The  current  function  fi(x,y)  of  the  flow  pattern  is  given  by 


After  transformation  to  the  coniooidal  coordinates  (X,fi)  this  is  easily 
integrated  to  give 

20(x,y)  ■ U - (■*  slnh  2p  •»  j|)  . (2.2) 

The  cavity  surfaoe  fil*,y)  =0  is  thus  tha  prolate  spheroid  p • 

./here 
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■ llCT*.)  • sinh  2u0  - Ho  , (2.3) 

and 


Ho  ■ m* . (2.4) 

Tho  axial  component  of  velocity  ia  given  by  the  aquation 

(2.5) 

In  tho  plane  of  ayanetry,  where  X * 0 it  follows  from  (l.l)  that  the  operator 

(yd/dy)  ia  equivalent  to  the  operator  ( -tanh  m d/d  jj  ) , ao  that  the 
suhatitution  of  (2.2)  in  (2.5)  yields  eventually 

(m-l)U«M“  tanh  u , (X  ■ 0).  (2.6) 

Hence,  using  the  relation  (1.6),  we  obtain  for  the  cavitation  nuafrer  the 
expression 

Q „G (no)  . [“-***  **  8^‘  ,.jj*  - 1 .(2.7) 

Lsinh  mo  oosh  uo  ~ u* I 

The  dimensions  of  the  cavity  are  obtained  immediately  from  (1.3),  thus 

y«/f  ■ ooseoh  Mo 

and  _ _ J (2.8) 

Xo/t  * ooth  Mo  * 

Moreover,  on  substituting  (2.6)  and  (2.7)  into  (1.8)  and  evaluating  the 
integral,  we  find  that 

■ slab*  Mo  (In  oosh*  Ho  - tanh*  uP).  (2.9) 


3.  Cavity  in  a stream  of  circular  croas-aootion 

In  the  previous  section  we  examined  the  flow  patterns  produced  by  a linear 
distribution  of  souroes  and  sinks  immersed  in  an  otherwise  uniform  unbounded 
stream.  The  object  of  this  section  will  be  to  examine  the  flow  patterns  pro- 
duced by  the  same  line  source -sink  imnersed  in  a stream  which  has  a finite 
circular  oross— section.  It  is  convenient  to  take  the  radius  of  the  cross- 
section  of  the  undisturbed  stream  as  unity,  so  that  certain  boundary  conditions 
have  to  be  satisfied  on  the  cylinder  y « 1 . 

He  «H«n  consider  two  cases,  firstly  when  the  stream  is  enclosed  by  a 
rigid  cylindrical  wall,  and  secondly  wham  it  is  opai  Ob  all  sides  to  an  atmosphere 
at  oonstant  pressure.  The  first  case  obviously  i closes  the  boundary  condition 
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that  the  cylinder  y « 1 shall  be  a stream  surface,  in  other  words  Q(x,1  )■ 
constant.  Moreover,  consideration  of  conditions  far  upstream  or  downstream 
3*  the  oavity  shows  that  the  constant  is  simply  \ U. 

In  the  second  oase  the  mathematical  boundary  condition  only  approximates 
to  the  physioal  boundary  conditions.  Par,  instead  of  ensuring  that  the  velocity 
magnitude  shall  be  constant  on  the  disturbed  free  surface , we  aotually  stipulate 
that  the  axial  ooqponant  of  velocity  ahall  be  oonstant  on  the  cylinder  y ■ 1 . 
The  validity  of  this  approximation  has  been  investigated  in  Appendix  I,  where 
it  is  shown  that  the  errors  involved  are  not  serious  in  that  part  of  the  field 
whioh  forms  the  subjeot  of  this  report, 

lbe  current  function  due  to  a unit  point  souroe  situated  at  (g,0) 
inside  s free  or  fixed  cylindrical  boundary,  having  the  equation  y ■ 1 , is 
given  by  (ref,  2,3  ) 

-fv  ♦ A(x-e,y,v,-0  when  x > ? > 

l (3.D 

or  jV  - - fl(S-*,y,v,-d)  when  x < £.  J 


Here 

- Jt  (*»,»)*  0xp^ ■*«/,«*)  . (3.2) 

k„  B is  the  m th.  positive  aero  of  J„(s),  and  v • 0 for  a free  boundary 
add  v ■ 1 for  a fixed  boundary. 

On  integrating  from  ( . -f  to  £ = f , to  obtain  the  current  function 
t^(x,y)  due  to  the  line  source-sink,  the  result  takes  different  analytical 
farms  aooordii^  to  the  relative  magnitudes  of  x and  f . The  two  essentially 
different  oases  are  as  follows: 

(i) 

*hen  0 < x < f , 

- f {-•}■*♦  0(x-C^r,v,-l)!  ed$*  f{f«'-y  ‘5-o(5-x,y,v,-i)|  E, 

1 -t  Jx 

■ [C°(*-C^r,v,-2)  ♦ n(x-5jr,w,-3)]if  ♦ f v(t * - x#) 


♦ (?0(5-«^,w,-2)  -0(e-xjt,w,-3)  ♦ f . (3.3) 


(ii) 

Whan  f < x 

ft 

l^(x,y)  ■ J ^ i-fv  ♦ ttfx-e»Jrfy,-«)l  K <*€ 

a (5  * n(x-jf^r,v,-3)]^f  (3.4) 

The  current  function  of  the  contained  flow  pattern,  0(r ,y) , is  simply 
the  sub  of  the  current  function  due  to  the  line  saarce-sink,  {^(x,y), 
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and  that  dua  to  tba  uniform  stream:  thus 


°(*,y)  ■ t U ♦ 0,(x,y)  . (3.5) 

Tha  function  0«(x,y)  is  tha  sum  of  an  infinite  series  of  products 
of  Bessel  functions  and  exponential  functions,  but,  by  means  of  certain 
transformations  listed  in  Appendix  II  , it  can  be  expressed  in  another  form. 
Thus, 


2(^(x,y)  « - sinh  2„  - „)  ♦ *■  0„  , f*  ♦ » (*,y.O  , (3.0 

where  GU  a are  constants  and  R (x,y,f)  ia  the  sum  of  an  infinite  series 
whose  nth.  tern  is  a homogeneous  polynomial  in  x,y  and  f of  degree  (n  + 3), 
This  farm  demonstrates  the  similarity  between  the  current  function  for 
unbounded  flow  (2.2)  and  that  for  bounded  flow  (3.5),  and  will  be  particularly 
usefhl  in  a later  section  for  developing  a first  order  theory  for  small 
oawities. 

The  equation  of  the  dosed  stream  surface  nay  now  be  obtained  by  equating 
D(x  jr)  to  sero.  The  resulting  functional  relation  between  x and  y is  rather 
complex,  but  we  shall  investigate  only  the  overall  dimensions  of  the  cavity, 
namely  x,  and  y0  . These  quantities  - the  semi-length  and  equatorial  radius 
of  the  oavlty,  respectively  - are  given  by  the  equations 


n(x„0)  - 0 andn(0,yo)  - 0 . (3.7) 


Next,  it  is  necessary  to  differentiate  y*fi(x,y)  partially  with 
respect  to  y , in  order  to  obtain  the  axial  velocity  conponent  nU  . <7e 

observe  that 


- 7 (V*tX*,y,«',p)|  ■ *(x,y,v,i»i)  , 


where 


d(*^r,»,P)  * J (-*„, »)P  «xp(”*iV,x)  . 

-iAKv»e' 

It  follows  from  (3.5)  that 

* y jy  [V*  n(x^r)^  - u - *B(x,y)  , 


(3.8) 


(3.9) 


(3.10) 


where 


#B(xjr)  ■ #(x-^,w,-l)  ♦ 

♦ [5  *(e-*,y,i',-1 ) - ^(5-x,y,v,-2)]f 

or 

* [5  ) ♦ *(fc-£,y,«',-2)]ff 


v(f*-x») 

(Oexcf) 


► (3.1i) 


If  < x)  J 
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On  taking  use  of  the  transformations  given  in  Appendix  II  , one  oen 
express  0*10)  in  the  fora 

(■-O  & - « - + cosh  X sinh  y { seoh(x+s)  ♦ seoh(X-w)j  ♦ * 0„  , f * ♦ R (x,y,f)} 

(3.12) 

alternatively , this  expression  may  he  obtained  directly  from  (3.6). 


In  particular,  at  the  equatorial  plane  x and  X vanish,  so  that,  frc 

(3.10,11,  12,) 

(n-1)  0 - 2f  #(fjr,»,,-l)  -*(f,j,„,-2)  ♦ 2*( 0jr,w,-2)  ♦ „ f». 


• u - tanhy  ♦ f C9ft  f»  ♦ R (0,y,f)  , 


> (XwO)  (3.13) 


and  the  corresponding  f omuls  for  m 0 is  obtained  by  simply  adding  the 
suffix  sero  to  the  symbols  m,  y and  y . Finally,  Q follows  from  equation 
(1.6),  and  (1.6)  then  yields  Cfo  after  a numerical  integration. 

4.  First  order  theory  for  tap-dimensional  oavities 

Before  developing  a first  order  theory  for  axially  symmetric  cavities 
it  seems  advisable  to  set  on  record  the  corresponding  results  for  two- 
dimensional  oavities  behind  a flat  plate.  The  exact  mathematical  theory  of 
a two-dimensional  cavity  behind  a flat  plate  in  a bounded  open  or  closed 
stream  has  been  developed  by  Simons  ( ref.  1 ) , but  a certain  amount  of 
algebraic  manipulation  is  necessary  in  order  to  extract  the  results  which 
we  require. 

(»)  Closed  tunnel 

A suitable  starting  point  far  developing  a first  order  theory  in  the 
case  of  the  closed  tunnel  is  furnished  by  Simmon's  equations  (47,  54,  56 
and  62)  (ref.  1)  . In  the  present  notation,  these  becomes 

T km^rc  * g(  6)  «•  k'  sec  6 tan ’ ’ ( k ' tan 5)  f ' 

•y  tanowxo  * k*  K sin  6 - sec  ^i-k'sin*®)  ,/#  |K  E(  6)  - B F(6)J  , 

4.kmo«yo  ■ g(6)  + k'  6 sec  6 , ^ 

kwcCi/lunp*  g(6)  , 

J 

where 

g(6)  » sec  6(l-k*sin*6),/t  Jk'  B(6)  ♦ (E'-*')  F(6)|  - k*  K'sjni,  (4.2) 


and 


ton0  « sec  6 {(l-4c*sin,6),/*  > k'j  , (4.3) 


■10. 


The  incomplete  elliptio  integral*  ere  defined  by 
B(6)  ■ j (l  -kSinSi)’ /f  du  ' 

► 

,6 

*6)  « f (l-k*sin«u)-,/,*i  , 

,•  e 

end  the  complete  elliptio  integrals  B end  K by 


(*►•4) 


B - B(*/2)  , K - l(#/2)  . (4.5) 


Primsd  letter*  for  oomplete  elliptio  integral*  indicate  that  the  modulus 
ha*  been  replaced  by  the  oo-oodulua  k '(kT+  k**  ■ l). 

The  assumption  that  k ' is  snail  leads  to  aerie*  expansions  valid  far 
small  value*  of  the  oavitation  number.  Thus,  from  (4.2), 


g(6)  ■ ^ rk"  tea  5 *eo  & jl  ♦ 0(k#,)J  , 
and  hano*  equations  (4.1)  take  the  farm 


•J-  km«eo  ««  ±(/r+K)  k'*  tan 6 *eo&  (l  ♦ 0(k'*)}  , ■> 

T a gd”6{l  ♦ 0(k'*) j , 

- 

4 kaovyo  ■ k'  6 eec5  Jl  ♦ ■£  k'  tan{/6  ♦ 0(k'*)]  , 
kwoCf/Umo  » k'*  tand  *eo6  jl  ♦ C^k'*)}  , 


and  (4.3)  beooaes 


km,  - 1 ♦ k'  s*o6  ♦ 0(k'*)  . (4.7) 


On  dividing  th*  last  three  of  aquations  (4.6)  by  the  first,  the 

of  the  oavity  are  axpraased  in  tarms  of  the  width  of  the  aavitating 
plat*)  thus 


S-W-TW! 

h ♦ 0(k'*)!  . 


(4.8) 
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The  equations  (4*8)  involve  the  two  parameters  ft  , which  broadly 
represents  the  effect  of  the  tunnel  sise,  and  k'  , which  broadly  represents 
the  effect  of  cavitation  number.  We  shall  gpw  separate  these  two  effects, 
so  far  as  this  is  possible,  by  writing  Xg,  y.  and  Op  far  the  limiting 
forms  of  xs,  yo  and  Cq  , respectively,  as  o tends  to  sero,  that  is  to  say, 
as  the  tunnel  cross-section  increases  indefinitely,  First  making  the 
substitution,  obtained  from  (4.7) , 


k'  seoft  wjQjl  - ^ Q ♦ 0(Q*)  j , (4.9) 


we  find  that 


is  _ .16  1. 1 Q ♦ 

c WS4  Q* 

Z*  - - _B  1 ♦ T Q ♦ 0(0*1 

o «r»4  w*4  Q 


(4.10) 


and 


S,  * Jl  ♦ Q ♦ 0(Qb)I  , 


where  the  right-hand  members  are  functions  of  Q only.  The  corresponding 
correction  factors,  expressing  the  effect  of  tunnel  size,  are  then  found  to 
be 


and 


i$sr  l’  * . 


y° 

Jo 


1SX  f1  ♦ 


♦ o(Q*)  . 


(4.11) 


It  will  be  noticed  that  these  correction  factors  are  not  entirely  independent 
of  the  cavitation  number,  although  when  the  cavitation  number  is  small  the 
inter-dependence  is  only  slight. 

The  precise  physical  significance  of  6 is  apparent  from  (4.6),  which 
shows  that  it  is  related  to  the  length  of  the  cavity  by  the  equation 

•^w~x0  * gd"1 6 |l  - £ Q ♦ 0(Q*)J  . (4.12) 
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Accordingly  we  plot  in  Fig.  3 the  behaviour  of  xJ%o  and  %/y9  aa 
function*  of  Xo  in  the  limiting  oase  when  Q ■>  0 , (o  is  then  negligible 
in  comparison  with  y0  or  y0>)  using  6 as  a parameter. 

Up  to  this  point  no  restrictions  have  been  imposed  on  the  parameter  J, 
which  is  free  to  range  from  zero  to  <1  w • Aooordingly  our  results  are 
valid,  for  small  values  of  Q , whatever  may  be  the  ratio  of  the  length 
of  the  cavity  to  the  width  of  the  tunnel.  It  should  be  noticed  however, 
that,  as  Simmons  pointed  out,  the  ratio  of  the  width  of  the  cavity  to 
the  width  of  the  tunnel  must  necessarily  be  small  when  Q is  small. 


tie  now  particularize  our  results  to  the  case  of  cavities  which  are 
short  compared  with  the  width  of  the  tunnel.  Then  & is  small  and  the 
equations  (4.11)  reduce,  by  virtue  of  (4.12),  in  the  oase  of  vanishingly 
small  cavitation  number,  to 


These  results  may  be  sinply  expressed  in  words  as  follows: 

(i)  For  small  values  of  the  cavitation  number  the  ratio  of  the  length 
of  a cavity  to  the  width  of  the  stream  may  be  of  any  order  of  magnitude, 
although  the  width  of  the  bubble  is  necessarily  small  compared  with  either 
of  these  dimensions. 


(ii)  The  true  drag  coefficient  of  the  cavity  tends  to  be  independent  of 
the  blockage  ratio,  for  snail  cavitation  numbers. 

(iii)  When  the  ratio  of  the  length  of  the  cavity  to  the  width  of  the 
stream  is  small,  and  the  oavitation  number  is  also  small,  the  length  and 
breadth  of  the  cavity  are  both  increased  by  an  amount  proportional  to  the 
square  of  this  ratio,  and  the  proportionate  effect  on  the  length  is  twice 
as  great  as  that  on  the  breadth. 


(b)  Free  Jet 


In  the  case  of  the  free  jet,  one  of  the  present  authors  has  adapted 
Simmons's  results  to  a form  which  furnishes  a suitable  starting  point  for 
the  present  treatment.  Thus,  from  equations  (3.8),  (3.9)  end  (3.10)  of 
ref.  4 we  have 


•fwBeo  ■ 4 q,/*  * q*  Ci(w')  , 

Jr  mtXo  ■4q,/*  * q*  x»(«')  , 

i *0 

f mstfo  ■ 4 q"*  2 q*  Yt(w') 

i «0 

and  <|>  troCj/mo  ■ 4 q,/#  2#  q*  Dt(  *')  , 


(4.14) 


where 

*'  * ♦ «)  (4.15) 

and  C»(«r')«  Xt(r'),  Y»(«'),  D*(ir')  are  known  functions.  In  on* 

respect  these  initial  equations  are  simpler  than  those  (4. 1,2,3)  for  the 
dosed  tunnel,  since  one  of  the  mathematical  parameters,  namely  *' , 
depends  on  the  cavitation  number  Q only.  On  the  other  hand,  of  course, 
they  possess  the  compensating  disadvantage  of  involving  infinite  series. 

The  other  mathematical  parameter,  namely  q,  represents  broadly  the  effect 
of  tunnel  aise  and  is  accordingly  analogous  to  6 in  the  previous  case. 

In  order  to  study  the  case  of  small  cavitation  numbers,  we  expand  each 
of  the  coefficients  of  q*  in  equations  (4.14)  as  series  is  ascending 
powers  of  k ' and  obtain  the  following  results 


and 


^wisoo  a (w  ♦ 4)  q,/f  T,(q)  *'*  |l  ♦ 0(  *'*)!  , 
JfimioXo  = 4 q’/#  Tg(q)  {l  «.  0(w'*)i  , 

•J-wm«yo  =>  q1'*  [4  r^q)  k'  ♦ v f,(q)  ♦ 0(*'*)] 

j-ircOD/mo  a 2 w q’/l  T,(q)  w'*  [l  ♦ 0(  *'*)}  , 


(4.16) 


where 


and 


r,(q)  = 1 ♦ 4 q * 6 q*  ♦ 8 q*  ♦ 0(q4)  , 
ra(q)  « i ♦ f-  \ ♦ $ q*  ♦ f q’  ♦ o(q4) 
r»(q)  * i ♦ f q ♦ f q’  ♦ £>(q4) , 


(4.17) 
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Division  of  the  latter  throe  equations  of  (4.16)  by  the  first  gives 
expressions  for  the  cavity  dimensions  in  terms  o'  the  width  of  the  oavltating 
plate,  thus 


and 


a 


h ♦ o(*'k)j . 


(4.18) 


The  similarity  between  equations  (4.18)  for  the  free  jet  and  (4.8)  for 
the  closed  tunnel  is  already  clear,  and  it  is  perhaps  worth  mentioning  in 
pacing  that  the  relation  between  k ' and  tc ' is 


*'  " t#nh  fa^T  t“h‘”  k'}  ’ 


(4.19) 


so  that  «'  and  k'  tend  to  equality  when  small. 

Proceeding  as  in  the  earlier  case,  we  now  write  x0,  y0  and  (5) 
for  the  limiting  forms  of  Xo  , y0  and  Op  , respectively,  as  q tends  to 
sero,  that  is  to  say,  as  the  tunnel  arose-section  increases  indefinitely. 
On  usii«  equation  (4.15),  this  leads  at  once _to  ajrestatement  of  equations 
(4.10),  shewing  that  our  two  definitions  of  x0,  yp  and  7^)  are, 
in  fact,  consistent.  The  corresponding  correction  factors,  expressing  the 
effect  of  tunnel  sise,  are  found  to  be 


and 


*a 

x0 


♦ 0(q»)j 


\ 


l1  * ««">1 


~ = 1 ♦ o(q*)  .•  J 

°D 


(4.20) 


See  note  at  end  of  this  section. 
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The  physical  significance  of  q may  be  deduced  from  the  second  of 
equations  (4. "16),  which  becomes 

•?  wxo  = 4 q1/2  tj(q)fl  - i Q + 0(Q2)|  . (4.2I) 


Using  q as  a parameter,  we  can  now  plot  in  Pig,  3 the  variation  ofjq,/^ 
and  y0/y0  as  functions  of  x0  , in  the  limiting  case  when  Q vanishes. 

The  above  results  are  true  for  quite  general  values  of  q between  0 and 
1 , although  coefficients  of  higher  powers  of  q than  are  contained  in 
(4.17)  would  be  necessary  to  obtain  reasonable  accuracy  in  the  evaluation 
of  the  functions  I\(q),  T2(q)  and  r3(q)  when  q becomes  appreciable  in 
size.  The  evaluation  of  further  coefficients  would  be  tedious  but  not 
difficult. 

If,  however,  we  now  particularize  to  the  case  of  small  q,  that  is 
to  say,  when  the  length  of  the  cavity  is  small  compared  with  the  width 
of  the  jet,  the  equations  (4.20)  reduce,  by  virtue  of  (4.21),  to  the 
formulae 


x</x0  = 1 - ~ xo2  + 0(x04).  ^ 


and 


yc/yo  = 1 - fg  x02  + 0(xc4)  f (q  - 0)  (4.22) 


W1 


These  particular  results  have  already  been  obtained  in  an  earlier  paper, 
but  the  general  results  (4.20,  21 ) are  new. 

It  will  be  observed  that  equations  (4,22)  are  very  similar  in  form 
to  equations  (4.1 3) > the  only  difference  being  that  the  coefficients  in 
the  case  of  the  free  jet  are  of  opposite  sign  and  of  half  the  magnitude 
of  those  pertaining  to  the  closed  tunnel, 

NOTE: 


It  will  be  observed  that  in  the  limiting  case  when  Q = 0 , equations 
(4,11)  and  (4.20)  both  imply  that  the  true  drag  coefficient  is  completely 
independent  of  cavity  size.  This  statement,  however,  has  not  a great  deal 
of  practical  significance  since  in  this  limiting  case  the  cavitating  head 
itself  must  vanish.  This  is  clear  from  equations  (4.6)  and  (4.16),  for  c 
vanishes  with  k'  (or  «'  ) , whatever  the  value  of  6 (or  q). 

At  first  sight  this  result,  as  regards  the  free  jet,  seems  to  be  at 
variance  with  previous  results  (see  for  example,  refs.  4 and  6),  where 
the  variation  of  the  true  drag  coefficient  at  vanishing  cavitation  number 
is  actually  studied  quantitatively.  There  is  in  fact,  however,  no  con- 
tradiction, since  in  the  earlier  papers  the  cavity  was  allowed  to  be  open 
at  the  rear,  whilst  in  the  present  paper  it  is  necessarily  closed. 
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5.  First-order  theory  for  axially  symmetric  cavities 


With  the  results  of  the  exact  two-dimensional  theory  to  act  as  a guide, 
we  are  now  in  a position  to  develop  the  first  order  formulae  from  the 
approximate  axially  symmetric  theory. 

Since  it  is  required  to  compare  the  properties  of  cavities  in  bounded 
and  unbounded  streams  at  a given  value  of  the  cavitation  number,  the 
fundamental  connection  is 


m0  = in  o 


(5.1) 


For  the  other  variables,  we  express  the  bounded  stream  value  as  the  sum  of 
the  unbounded  stream  value  and  a small  increment.  Thus,  for  example,  the 
equation  (3»13)  evaluated  at  the  equator  of  the  cavity  becomes 


(mo  - l)(U  + ty)  = do  + Ado  - tanh  do  - sech*  do  Ado  + §■  Cv  a f3  + R (5.2) 


where  R will  be  used  to  denote  a remainder  composed  of  smaller  order  terms. 
On  the  other  hand,  equation  (2.6)  evaluated  at  the  equator  of  the  cavity 

becomes 

(m0  -1)  U = do  - tanh  "do  , (5.3) 

and  accordingly,  by  subtraction, 

(*0  " 1 ) AU  = tanh2  do  Ado  + -f-  Cu> 2 f 3 + R . (5.4) 

Similarly,  since  the  value  of  d at  the  point  (o,yo)  is  do  > we  find 
from  equations  (3.5,  6 and  7)  that 

U + * £ sinh  2 do  + cosh  2 d0  Ado  - do  - ado  - f C , f3  + R , (5.5) 

whilst,  from  (2.3)> 

U = 4 Sinh  2 do  -do  . (5.6) 


A subtraction  now  yields  the  result 


AtJ  = 2 sinh*  do  ^do  - $ Cv  2 t 3 + R . 


(5.7) 
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Again,  the  value  of  ^ at  the  point  (xo,0)  is  y^i.e.  -f  , 

so  that  (3.5,  6 and  7),  together  with  (2.3,  4)  give  a result 
.recisely  siriular  to  (7.7/,  but  with  the  subscript  «•  replacing  tnc  subscript 
o throughout,  Since,  however,  we  have  already  noticed  (2.4)  that  ~o  = "73 „ 
it  follows  that  u0  = Ua>  » and  the  cavity  in  the  boun'.ed  stream  is  still 
a prolate  spheroid,  to  the  present  order  of  accuracy,  ibis  fact  is  also 
clear  directly  from  (3.6). 

Elimination  of  U between  equation  (5.3)  and  (5.6)  yields  an  expression 
for  (m0  - l)  , which  may  be  equated,  to  a similar  expression  obtained  by 
eliminating  AU  from  (5.4)  and  (5.7),  thus 


_ ho  - tanh  un  _ tanh2  Un  4/n+  J Cy.tf3  + R (5.8) 

f sinh  2 q0  - do  2 sinh2  Uo  Ado  - -j  Cy  2 f3  + R 

Now,  when  the  cavitation  number  is  very  small,  these  two  expressions  for 
(mo  - l)  must  tend  to  zero  and  hence  ~jl0  must  tend  to  infinity,  both  the 
numerator  and  the  denominator  of  the  left-hand  member  of  (5.8)  are  then 
clearly  dominated  by  the  leading  terms.  It  may  then  be  deduced  that,  in 
the  right-hand  member  of  (5.8),  the  numerator  is  lominated  by  the  second 
term,  whilst  the  denominator  is  dominated  by  the  first  term.  It  follows, 
since  sinh  2 Uo  ~ 2 sinh  2 Uo  , that 


2 ho  A^o  ~ Cyf2  f3  (Do  - . (5.9) 

In  order  to  deduce  the  corresponding  change  in  the  proportions  of 

the  cavity,  observe  that,  from  (1.3), 


y0/x0  = sech  ho  , 


so  that  _ _ tank  y 0 Aq0  , (5.10) 

y<A  0 


Combining  (5.9)  and  (5.10),  we  find  that 

~ _ £i^£!  (ho-»).  (5.11) 

yo/^o  3 Uo 

The  numerical  values  of  - Cv  2 are  given  in  Appendix  13  as 

- Co  2 = .205911 ...  -c,  2 =-  .79632)  ...  , (5-12) 

* t ' 

so  that  the  corrections  to  be  applied  in  a closed  tunnel  are  of  opposite 
sign,  and  roughly  four  times  in  magnitude,  those  to  be  applied  in  a free 
jet.  This  makes  a reasonable  comparison  with  the  two-dimensional  case, 
where  the  corresponding  factor  is  exactly  two. 
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Again,  the  corrections  to  be  applied  are  urcp crt ioiial  to  the  cube 
of  the.  ratio  of  the  length  of  the  cavity  to  the  width  of  the  axially 
symmetric  jet,  whilst  in  tv. „ dimensions  the  corrections  were  proportional 
to  the  square  of  this  ratio. 

There  is  an  important  dissimilarity  between  the  two-dimensional  and 
axially  symmetric  results,  hoyrever,  in  that  the  numerical  coefficient  in 
the  latter  case  tends  logarithmically  to  zero  with  the  cavitation  number, 
whereas  in  the  two-dimensional  case  it  tends  to  a finite  limit. 

6.  Use  of  the  drag  coefficients  to^determine  cavity  dimensions 

The  width  of  the  cavitating  plate,  namely  2c,  does  not  appear  in  the 
analysis  which  we  have  developed  for  the  axially  symmetric  cavity,  so  that 
it  is  impossible  to  determine  unaided  the  variation  of  the  actual  cavity 
dimensions,  as  distinct  from  the  thickness  ratio  yo/xo  • However,  it  is 
clear  from  the  two-dimensional  theory  that  the  true  drag  coefficient  Cp  , 
based  on  the  cavitating  plate  v/idth,  is  very  insensitive  to  boundary  effects, 
and  by  assuming  that  the  same  is  true  in  the  axially  symmetric  case  we  can 
calculate  the  variation  of  x0  and  y 0 separately. 

Prom  equations  (1.7,  8),  the  true  drag  coefficient  is  given  by 


C 


D 


(y02/c2) 


(6.1) 


•where 


= o - sinh2  do  l(  g0)  , 


(6.2) 


and 


l(do)  = 2 f 


(m  - 1 )z  y dy  . 
y o 


(6.3) 


Taking  logarithms  in  equation  (6.1)  and  subsequently  performing  the  operation 
A , we  find  that,  since  Cp  = Cp  by  assumption, 


2 -&fl  = . ^ 

% c r • 


(6.4) 


The  cavitation  number  Q is  the  same  in  both  bounded  and  unbounded  cases, 
so  that 


- ACp  = sinh  2 do  A^0  l(u0)  + sinh2  do  A I (do)  . (6.5) 


On  making  the  transformation  to  the  conicoidal  coordinates  (l.l)  and  remember- 
ing that  the  path  of  integration  is  along  X = 0 , (o.3)  becomes 

do 

l(do)  = 2 f (m  - l)2  cosech2  u cotli  n dd  , (6.6) 

•'sinh*1  f 
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whence,  in  particular, 


l(Mo)  = 2 Vo  coth  /j0  - 2 In  cosh  y0  - p0*  cosech1,  fj0  - 1,(6.?) 

and 

l(Fo  ) ~ (In  4 - 1 )/u  * (u o - «)  • (6.7a) 

The  quantity  Al(^0  ) may  be  evaluated  from  (6.6)  as  the  sum  of  four 

parts,  thus 

dl(  Vo  ) = d,  + A,  + As  + A4  f (6.8) 


where  d,  = 2 (m0  - l) * cosech*  ji0  coth  Afj0  + R 


r f * 
s c v,ar 

U* 


(m0  - 1 ) 


( Mo  ®)  (6.9) 


and 


d. 


f 

(m  - 1 ) * cosech*  (i  coth  //  d/j  + R 


( Vo  - ®)  « (6.10) 


The  two  contributions  A,  and  Aa  result  from  the  alteration  of  the  upper  and 
lower  limits  of  integration  in  l(  pe)  . The  remaining  contributions,  d3 
and  4* , result  from  applying  the  operator  A to  the  integrand.  For  a given 
value  of  (i  , the  only  factor  in  the  integrand  which  is  susceptible  to  A is 
(m  - l)*  . Now,  from  (3.13)  and  (2.6) 


f4 

1®* 


df(m-1  ) U|  = U A(m-1  ) + (m-1  ) AU  = f Cv>2  f3  + R 


so  that  U A(m-1)  ~ CVfJ  fs  - (m-1  ) AU  . 


(6.11) 
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On  using  this  result  we  now  find  that 


cosech  2 u coth 


(j4)i  + R 


fc  v af  3 

~~r~ — (ho  - *>)  ; (6.12) 

u 2 


and 


A4  = ^2  l(„0)  + R 


f 0V  i f 3 1 

3 y*  ho" 


(ho  **  »)  . 


(6.13) 


It  is  clear  that  when  f is  small,  and  when  /j0  . «°  that  is  to  say, 

when  the  cavitation  number  is  also  small.  A,  A8  and  A4  are  all  small 
compared  with  A,  , so  that  equation  (6.5)  becomes 


(6.14) 


The  first  term  inside  the  square  bracket  may  be  neglected  by  comparison 
with  the  second  term,  when  ~Ji0  is  large. 

Examining  the  relative  magnitude  of  the  terms  in  (6.2),  we  find  from 
(5.8)  that 


Q ~ 2(m,-1 ) ~ 2 h</sinh2  ho  (ho  *«>),(6*15) 

and  from  (2,3)  and  (6.7a)  that 

sinh2  ho  l(ho)  ~ (In  4 - 1_)^osh2  "ho  (ho  * •>);(^*1^ 

so  that  the  second  item  is  negligible  compared  with  the  first. 

Substitution  of  the  results  (6.14>  15  and  16)  in  equation  (6.4) 
produces  the  result 

Ayo  fcy.2f3  sinh»70 
To  ~ coslf -j;0  2 7T0 


cy,af3 
3 ho 


(4o  - •)  . (6.17) 
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Finally,  since 


In  x'o  = In  y0  - In  (yc/x0)  , 


we  find,  on  applying  the  operator  A , that 

*£&  = J&L&  Kyo/xp) 

xo  yo  yi/^o 


2CV  tf3 

— ■ lUk 

3 Wo 


( Wo  •*  “),  (6.18) 


by  virtue  of  (5.11)  and  (6.17). 

The  results  (6.17)  and  (6.18)  for  the  axially  symmetric  cavity  are 
analogous  with  those  (4.13  and  22)  for  the  two-dimensional  cavity.  They 
state,  in  fact,  that,  for  a small  axially  symmetric  cavity  at  low  cavita- 
tion number  in  a closed  or  open  tunnel  of  circular  cross  section,  the 
boundary  effect  is  proportional  to  the  cube  of  the  ratio  of  cavity  length 
to  tunnel  breadth,  and  the  proporticnate  effect  on  the  cavity,  length  is 
twice  as  great  as  that  on  the  cavity  breadth.  Furthermore,  the  boundary 
effect  is  a shrinkage  in  the  case  of  the  open  jet  and  an  expansion  in  the 
case  of  the  closed  tunnel,  the  latter  effect  being  roughly  four  times  as 
great  as  the  former. 


7.  Numerical  results 


The  simple  first  order  results  deduced  in  sections  5 and  6 apply  only 
to  the  doubly  limiting  case  when  both  Q and  f (or  x0  ) are  very  small.  One 
of  these  limitations,  that  of  small  cavitation  number,  applies  also  to  the 
more  general  theory  developed  in  sections  2 and  3,  for  the  prolate  spheroids 
investigated  therein  satisfy  the  constant  pressure  condition  of  cavitating 
flow  only  when  the  thickness  ratio  y0/x0  is  very  small. 

The  condition  that  the  length  of  the  cavity  should  be  small  compared 
with  the  width  of  the  tunnel,  however,  is  not  necessary  to  ensure  the 
validity  of  the  theory  of  sections  2 and  J>.  It  is  true  that  the  width  of 
the  cavity  must  be  small,  in  order  that  the  mathematical  boundary  condi- 
tion applied  to  the  surface  of  the  free  jet  shall  so.tisfactorily  represent 
the  actual  physical  boundary  condition.  But  for  any  given  cavity  length, 
the  width  of  the  cavity  is  bound  to  be  small  when  the  cavitation  number 
is  small  enough,  so  that  we  do  not  need  to  stipulate  that,  in  addition, 
the  length  of  the  cavity  shall  be  small. 

Accordingly,  it  seems  worth  while  to  evaluate  the  results  of  the 
general  theory  for  a wide  range  of  values  of  the  non-dimensional  cavity 
length  x0  , and  for  various  small  values  of  Q.  In  order  to  be  convenient 
for  practical  use,  the  results  have  to  present,  for  given  values  of  Q, 
the  variation  of  the  thickness  ratio  yo/xo  as  a function  of  some  easily 
measured  length,  such  as  x0  • In  view  of  the  fact  that  Q is  only  obtained 
in  terms  of  x0,  y0  and  f after  a fairly  lengthy  calculation,  some  care 
is  necessary  in  order  to  reduce  to  a minimum  the  labour  and  inaccuracy 
involved  by  inverse  interpolation.  The  method  of  calculation  which  was 
eventually  selected  will  be  briefly  described. 
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The  oavity  in  a bounded  stream,  at  a given  value  of  the  oavitation 
number  Q,  la  aharaoteriaed  by  a trio  of  value  a jis,  y*  and  f , oorv- 
neoted  by  the  relation  (1.3).  Thia  value  of  oorresponds,  in  an 
unbounded  atraaa,  to  a slightly  different  value  of  the  cavitation  number, 
which  we  ahall  denote  by  Qw  * G(«e). 

Similarly , the  oavity  in  the  bounded  atrean  ia  aharaoteriaed,  for 
given  Q,  by  the  triQ.ji«.»  x,  and  f,  oonneoted  by  the  relation  (l.j),  #• 
aooordingly  define  Qx  ■ G(y  J . , that  ia  to  aay,  the  oavitation  matoer  in 

the  unbounded  stream  whioh  would  yield  a oavity  of  the  length  aotually 
observed  in  the  bounded  stream. 

The  essence  of  our  method  ia  to  form  estimates  of  Qy  and  Qx  frcai 
the  given  value  of  Q and  henoe  by  using  (2.7)  to  calculate  pa  and  , 
which  yield  the  required  values  of  y,  and  x..  The  process  is  iterative, 
since  each  estisnte  of  y#  and  x,  enables  a better  eatlsation  of  Qs 
and  Qx  to  be  made,  and  these  in  turn  lead  to  improved  estimates  of  y, 
and  x*. 

The  first  step  was  to  tabulate  Q as  a function  of  Tie  > using  (2.7), 
at  intervals  sufficiently  close  to  oaks  inverse  Interpolation  simple  and 
accurate  (Table  I).  Next,  taking  in  turn  the  values  y » 0.05  (0.05)  0.35, 
combined  with  the  values  f - 0.5  (0.5)  2.0,  the  functii  • U(y0,  f)  was 
calculated  from  (3.5)  by  setting  0(0, y,)  - 0 . For  the  same  values  of 

y8  end  f,  the  product  (a,  - l)  U was  also  calculated  from  eqiation 
(3.13).  Sisf>le  division  then  yielded  m„  , ftrom  which  Q followed  by  virtue 
of  (1.6).  Far  the  seme  ratios  yjt  , the  function  Q„  ^s  calculated 
directly  from  (2.7).  fith  f as  a parameter,  the  difference  Q„  - Q was 
tabulated  as  a function  of  y. , and  the  resulting  second  differences  were 
sufficiently  steady  to  ensure  satisfactory  direct  interpolation  (Table  11). 

In  a similar  manner,  but  setting  Q(ao,0)  ■ 0 in  place  of  0(0,y*)  ■ 0, 
the  function  Qx  - Q was  tabulated  as  a function  of  x*  , "1th  f as  a para- 
meter, intervals  of  0.02  in  the  argument  being  quite  small  enough  to  ensure 
accurate  interpolation  (Table  III). 

Par  a riven  valug,  of  Q and  one  of  the  tabular  values  of  f , the  first 
estimate  of  7<j  (or  Q x ) was  simply  Q itself.  The  corresponding  value 
of  ne  (or  Umm  ).  obtained  by  inverse  interpolation  from  Table  I,  yielded 
a first  estimate  of  y0  (or  x0  ).  The  corresponding  estimate  of  the 
difference Tu  - Q (or  Qx  - S ) was  interpolated  from  Table  II  (or  III) 
and  an  ixqproved  estimate  of  Q ^ (or  Qx  ) thus  obtained.  The  cycle  was 
repeated  until  the  value  of  y 0 (or  x,  ) so  obtained  ms  not  subject  to 
further  change,  whin  J.t  was  assumed  to  be  correct.  Convergence  of  the 
process  was  rapid. 

The  advantage  of  this  method  of  computation  is  that  the  use  of  the 
unwieldy  formulas  involving  infinite  series,  which  ere  appliaeble  to 

streams,  la  kept  to  an  absolute  minimum.  The  only  inverse  inter- 
polation required  ia  for  the  ocsparwtively  simple  relation  (2.7),  *ich 
applies  to  an  unbounded  stream. 

The  calculations  of  Of)  were  planned  on  similar  principles,  this 
being  rendered  even  more  i^ortant  owing  to  the  faot  that  numerical  inte- 
gration was  involved.  The  expression  (1.8)  for  0f>  **  split  into  parts, 
thus 


$ 


(7.1) 


Og  ■ Q - (A  - B)A* 


where 


A a 


2 TV 


f (■  - 1 )*  If*  y 
*»« 


> 


an*  B ■ 2 yj* 


- (m 


O*  &“i  r . j 


(7.2) 


Dm  wall  numerical  Integral  B w tabulated  (Table  7)  for  the  saaa 
values  of  y#  sal  f as  those  used  in  the  compilation  of  Table  II. 
Then  for  one  of  the  ohcaen  values  of  Q,  and  a tabular  value  of  f,  the 
previously  oosgaited  value  ofy*  enabled  the  value  of  B to  be  inter- 
polated, whilst  the  value  of  A was  readily  oaloulated  analytically. 


The  results  of  our  computations  are  displayed  in  Pigs.  4-6. 

Pig.  4 shoes  tbs  behaviour  of  ths  ratio  (yB/x,)/(y^xB)  as  a function 
of  x#  , in  ths  aaae  of  the  open  tunnel  or  circular  cross- seotion.  A 
corresponding  curve  for  the  tso-diaanaional  ease  is  shown  dotted  for 
ooaparlson.  Pig.  5 shows  the  sans  function  plotted  against  the  blockage 
ratio,  that  la  to  say,  y,  in  the  oase  of  the  tao-dlnenaional  curve,  and 
y*  in  the  axially  ajmMtrio  caae.  Pig.  6 shows  the  effect  of  the 
Ares  stress  boundary  on  the  nodi  fled  drag  coefficient  baaed  on  the  oavity 
disaster.  As  explained  in  Seotion  6,  ee  aay  draw  the  Inference  that  the 
sextain  cross-sectional  area  of  the  oavity  behaves  in  an  approximately 
reciprocal  wanner. 
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•hen  this  report  was  in  course  of  preparation,  it  was  learnt  that 
a closely  siadlar  investigation  was  being  conducted  at  another  Goveraeent 
establishment,  neatly,  the  Admiralty  Besearoh  Laboratory,  Tedding  ton. 

It  was  therefore  agreed  to  divide  the  field  for  nunerioal  work  between 
the  two  establishments.  Accordingly,  the  numerical  results  depicted 
in  the  graphs  at  the  end  of  this  report  refer  only  to  open  Jet  tunnels, 
which  are  of  more  particular  interest  to  this  establishment.  It  is 
understood  that  a report  containing  similar  results  for  closed  tunnels 
is  being  published  by  A.B.L. 

8.  Conclusions 

The  mall  oorreotloaa  neoessary  to  measurements  of  axlslly  symmetric 
cavities  in  circular  tunnels  of  free  Jet  type  have  been  displayed  graph- 
ically. These  results  are  applicable  to  aasdl  values  of  the  cavitation 
number  but  to  arbitrary  values  of  the  ratio  (length  of  oavity) /(width 
of  aross-seotioo  of  tunnel)  .-x,  . 

Par  the  oase  tfien  x,  la  smell  ea  well  as  the  cavitation  number 
have  demonstrated  ths  following  results* 
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(a)  In  two-dlnenaiooal  flow 

(i)  the  fractional  inareaaes  (dosed  tunnel)  or  deoreasea  (Area 
Jet)  in  the  length  and  width  of  the  oarity  are  both  proportional 
to  i,',  The  factors  of  proportional! tar  are 

length  width 

oloaed  tunnel  w*/l2 

fr ee  Jet  - w*/24  . w*/48 

(ii)  the  true  dreg  ooeffioient,  baaed  on  the  oroee-eeotlonal  area 
of  the  oawitating  head,  is  unaffected  by  the  tunnel  boundary. 

(b)  X»  axjiynetrio  flow 

(i)  um  fractional  deorease  (oloaed  tunnel)  or  inareaae  (free 
Jet)  in  the  thickness  ratio  of  the  oarity  is  proportional  to  Jt,’* 

The  faotars  of  proportionality  tend  to  aero  with  the  oa vltation 
nusber,  but  they  are  in  the  following  ratio: 

oloaed  tunnel  - .796824 

free  Jet  ♦ .205911 

(ii)  assuning  that  the  true  drag  coefficient  ia  unaffected  by  the 
boundary,  then  the  fractional  increases  (closed  tunnel)  or  decreases 
(free  Jet)  in  the  length  and  width  of  the  oarity  are  both  proportional 
to  x ,*  . The  factors  of  proportionality  tend  to  sero  with  the 
cavitation  number , but  they  are  in  the  ratio: 


length  width 

closed  tunnel  2 x .796824  .796824 

free  jet  - 2 x .205911  - .205911 
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Appendix  I 


The  validity  of  tha  approximate  boundary  oondltioo  at  the  free  surface 

The  physioal  boundary  oonditlon  at  the  free  aurfeoe  of  the  Jet  differs 
frca  the  approximate  mathematical  boundary  oonditlon  In  two  respects,  tfe 
have,  In  faot,  aatiafied  the  wrong  boundary  condition  at  the  wrong  place, 
and  it  is  neeeaeary  to  see  whether  the  inaccuracies  involved  are  serloua. 

Owing  to  the  exoe  salve  amount  of  numerical  work  which  would  be  involved 
in  an  exaot  investigation,  we  shall  restrict  ourselves  to  a numerical  check, 
in  the  plane  of  symmetry , together  with  a survey  of  the  order  of  magnitude 
only  eleewhere. 

One  of  the  inaccuracies  involved  is  due  to  the  displacement  of  the 
free  boundary,  and  this  is  obviously  a maximum  at  the  plane  of  symmetry, 
where  we  ehall  assume  that  the  radiue  of  cross-section  is  i ♦ < . 

The  other  inaccuracy  is  due  to  assuming  that  the  velocity  magnitude 
at  the  boundary  is  adequately  represented  by  the  axial  oonponent  of  the 
velocity.  This  particular  error,  whiah  vanishes  at  the  plane  of  symmetry, 
is  obviously  a "cosine  effect"  of  the  angle  of  inoliSatlcn  which  the  free 
surface  makes  with  the  axle.  The  order  of  magnitude  of  the  maximum  angle 
involved  may  be  taken  as  given  by  e /f . 

The  equation  of  the  "free  boundary"  of  the  jet,  when  calculated  by 
the  present  approximate  msthbd,  is 


so  that 


* V - - ns(0,l)  Jl  ♦ 0(e))  . (i.O 

Neglecting  the  factor  jl  + 0(*)}  in  the  right-hand  member  of  (i.l)  we 
have  calculated  the  displacement  « as  s function  of  cavitation  number 
and  cavity  length,  the  results  being  shown  In  Fig.  7. 

A more  exact  theory  than  the  one  used  in  this  report  mould  have  to 
apply  the  free  boundary  oonditlon  at  the  free  surface  itself,  that  is  to 
say,  on  a surfaoe  lying  between  y * 1 and  y = 1 + « . It  may  be 

seen  Intuitively  that  the  results  of  such  an  improved  theory  for  a given 
value  off,  would  lie  somewhere  between  our  present  results  for  that  value 
off  and  the  corresponding  results  when  f is  replaced  by  f/fl  ♦ e) 

In  other  words,  the  plotted  points  in  the  Figs,  k - 6 are  a little  tc  the 
left  of  their  correct  position.  The  horisontal  coordinates  should  be 
increased  by  a factor  which  lies  between  1 and  (l  ♦ «) , [( 1 + Zt)  in 
the  case  of  Fig.  5j.  The  smell  values  assumed  by  e in  Fig.  7 make  it  clear 
that  the  correction  involved  would  be  small. 


The  boundary  condition  which  has  actually  been  applied  in  this 
report  ensures  that  the  axial  component  of  velocity  shall  be  constant  on 
the  cylinder  y = 1.  Owing  to  the  displacement  of  the  Jet  boundary, 
however,  it  is  clear  that  this  axial  coEg.  or.f.it  of  velocity  ./ill  not  be 
constant  on  the  jet  boundary.  Accordingly  v/e  shall  investigate  the 
difference  between  the  velocity  at  infinity  and  the  axial  co^w—.t  of 
velocity  on  the  equator  of  the  jet  boundary . Ideally,  of  course,  this 
cliff erence  should  be  zero,  since  the  axial  component  is  at  tliis  point 
equal  numei'ically  to  the  velocity  magnitude , 

From  equation  '3.10)  the  actual  difference  is  fives  by 


(1  - m)  0 = e (0,1  ) jl  + 0(e)]  . 


Now,  from  the  definitions  (3.2)  and  (3.9),  it  con  be  shown  that 
^ (x,y,v,p)  = y 0(x,y,r,p+l)  , 

from  which  it  follov/s  that,  watt  inf  v - 0 for  the  fr^e  boundary, 

* ^ CO.1'  = - f n(f.' A°)  + n(f,1,0,-l)  - 0(0,1  ,0,-1  ) . (1.3) 

Substituting  this  formula  in  (1.2)  and  using  the  previously  calculated 
values  of  e,  v/e  can  now  plot  (l  - m)tT  as  a function  of  Q and  x0,  the 
results  being  shown  in  Fif.  8. 

In  order  to  provide  a.  basis  of  cotqarieon  v/e  al30  show  in  Fif.  8,  by 
a dashed  line,  the  additional  contribution  to  the  axial  compon :nt  of 
velocity  at  the  centz'e  of  the  jet  due  to  the  presence  of  the  free  surface. 
According  to  equations  ( 1.6)  ard  (j.13)  this  nay  be  ° dressed,  after  some 
.manipulation,  as 


(m  - 1 ) 


U 


(x.y)  = (0, 


n) 


In  2 f - 1 - 2f  #(f  ,0,0, -1  ) + 2 <f>{  f,0,0,-2) 

(IA) 


It  is  clear  from  Fig.  8 that  the  variation  of  the  axial  component 
of  velocity  on  the  free  boundary  is  of  a smaller  ci'der  than  the  variations 
produced  at  the  jet  centre  by  the  presence  of  the  free  boundary.  In  fact, 
the  two  effects  differ  by  a factor  of  at  leo3t  5 in  all  the  cases  used  in 
constructing  the  curves  for  Fig3.  b - 6,  The  centre  of  the  jet  is  chosen 
as  the  point  of  the  basis  of  conqiarison  because  the  effect  of  the  free 
boundary  is  least  there,  so  that  the  difference  in  order  of  the  two  effects 
v/ould  be  even  more  marked  elsewhere. 

Finally,  it  i3  necessary  to  examine  the  effect  of  the  angle  of  inclina- 
tion between  the  jet  boundary  and  the  axis  of  symmetry.  This  results,  of 
course,  in  a discrepancy  between  the  axial  component  and  the  magnitude  of 
the  velocity.  Since  this  is  ft  cosine  effect  the  oi-der  of  magnitude  is  given 
by  the  square  of  the  order  of  magnitude  of  the  angle  of  inclination, 
namely  (e/f)*.  This  is  compared  with  the  .uantity  (l  - m) , measuring  tile 
variations  in  axial  component,  in  fig.  9.  Clearly  the  two  effects  are  of 
the  same  order,  and,  since  the  latter  effect  has  been  shown  to  be  unimportant, 
so  is  the  former. 
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AFH5NDIX  II 


It  has  boon  shown  (ref.  5)  by  one  of  the  present  authors  that  infinite 
series  of  the  types  appearing  in  equations  (3.2)  and  (3.9)  may  be  trans- 
formed into  other  serids  which  converge  rapidly  near  the  origin.  The 
particular  results  required  in  the  computation  of  equations  (3.3)»  (3*4) 
and  (3.11)  are  listed  below. 


J.0 

* -7  |r(Us)j  ♦ 


%**•••*)  - t Mi*.))"  - „j. 


-itta  Wi..)!  -}(i-.V(n.)]  - 

where  r = x ♦ y , s = x/r 

and  Pn(a)  is  the  Legendre  polynomial  of  order  n. 

The  coefficients  C v „ have  the  following  numeric:.!  values: 


c.,-1 

= 0 

c.,-, 

» - .375000 

O 

o 

o 

= * .435345 

° 

- ♦ 1.106824 

CV 

= 0 

C 1,1 

= - 1 

c0,» 

= - .205911 

C ,,, 

= .796824 

C0,4 

* ♦ .658857 

C »,  4 

= - 1.200470 

°»»» 

= - 5.14657 

0 i , • 

= ♦ 7.45829 

°o,. 

= ♦ 73.7500 

c ’»• 

= - 96.2205 

0 e , i e 

= - 1632.77 

C i,  to 

* ♦ 2070.91 

0 (n  > 0). 
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